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Abstract:  A  microchannel  reactor  for  hydrogen 
peroxide  decomposition  is  being  developed  for 
integration  with  fuel  cell  systems  that  can  power 
undersea  vehicles.  However,  the  catalytic 
decomposition  of  H202  is  predisposed  to  thermal 
runaway.  A  micro-scale  packed  bed  reactor 
(MPBR),  theoretically  capable  of  inhibiting 
thermal  runaway,  is  under  development  in 
COMSOL  to  illustrate  thermal  management  and 
oxygen  production  during  this  reaction.  The 
COMSOL  model  solves  mass,  energy,  and 
momentum  balances  to  simulate  temperature  and 
concentration  profiles  within  the  reactor.  Using  a 
stainless  steel  block  around  the  capillary  to  act  as 
an  extended  surface  for  higher  heat  removal  rates 
and  an  initial  volumetric  flow  rate  of  2e"9  m3/s 
(about  0.1  ml/min),  a  temperature  rise  less  than 
8  K  was  simulated  and  an  outlet  concentration  of 
1716  mol/m3  oxygen  was  achieved.  Overall  the 
results  indicate  that  thermally-controlled  oxygen 
generation  from  hydrogen  peroxide 
decomposition  is  feasible  in  a  microreactor 
provided  there  is  sufficient  external  surface  area 
to  facilitate  convective  cooling. 

Keywords:  Micro-reactor,  Oxidant,  Multiphase 
Modeling. 

1.  Introduction 

A  major  challenge  facing  the  design  of  air- 
independent  power  systems  is  continuous  energy 
generation.  High  energy  density  and  potentially 
low  operating  costs  make  fuel  cell  power 
systems  an  attractive  option  for  high  endurance, 
air-independent,  undersea  vehicle  applications.1 
However  air-independent  fuel  cell  power 
systems  require  an  oxidant  supply.  Hydrogen 
peroxide  (H202)  decomposes  into  water  (H20) 
and  oxygen  (02)  [H202  —>  H20  +  V2  02] 
providing  a  dense  source  of  oxygen  per  unit 
volume  that  makes  it  a  valuable  commodity  as  an 
air-independent  fuel  cell  oxidant. 

The  high  oxygen  density  (0.0215  02  moles/ 
ml  in  50%  H202  liquid  solution  compared  to 
0.035  02  moles/ml  in  solid  NaC103  for  example)2 
ease  of  handling,  and  commercial  infrastructure 


of  H202  add  to  its  appeal  as  an  oxidant  source.3 
Despite  these  advantages,  it  is  well  established 
that  catalytic  H202  decomposition  is  susceptible 
to  thermal  runaway,  which  has  historically 
limited  its  application  as  a  power  system 
oxidant.4  Microchemical  systems  by  virtue  of 
microscale  geometry  (relevant  fluid  dimension  in 
subunits  <  1  mm),  possess  high  surface-to- 
volume  ratios  resulting  in  heat  and  mass  transfer 
coefficients  capable  of  inhibiting  thermal 
runaway.5  Thus,  microchemical  systems  present 
a  theoretical  mechanism  to  harness  the 
aforementioned  benefits  of  a  hydrogen  peroxide 
oxidant  source,  and  at  the  same  time  prevent 
thermal  runaway  during  the  exothermic  reaction. 

The  aim  of  this  on-going  modeling  effort  is 
to  demonstrate  oxygen  production  and  thermal 
management  feasibility  during  the  catalytic 
multiphase  decomposition  of  hydrogen  peroxide 
in  one  of  the  subunits  of  a  microchemical  reactor 
system.  The  basis  of  the  present  model  is  a 
microchannel  reactor.  The  model  description, 
governing  equations  and  assumptions,  and 
boundary  conditions  follow  in  section  two.  Next, 
the  simulation  results  illustrating  temperature 
and  concentration  profiles  are  shown.  Finally,  a 
summary  is  given  for  the  current  model’s 
capabilities  and  limitations  and  areas  for  future 
modification  intended  to  better  evaluate 
multiphase  microscale  flow  effects  on  this 
reaction  are  identified. 

2.  Methods 

2.1  Model  Description 

In  order  to  simulate  the  catalytic  multiphase 
decomposition  of  hydrogen  peroxide  in  a 
microreactor,  a  channel  of  sub-millimeter  radius 
(0.5  mm)  containing  catalyst  was  modeled 
(Figure  1).  The  cross  sectional  geometry  of  the 
reactor  channel  resembled  a  half  moon.  The 
channel  proceeded  straight  down  the  length  of 
the  reactor  (5  cm).  Surrounding  the  channel  was  a 
rectangular  stainless  steel  316  block  (7x2x0.4 
cm).  To  mimic  upcoming  experimental 
conditions,  a  Plexiglas  (PMMA)  cover  (7x2x0. 1 
cm)  served  as  the  microreactor  seal. 
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Figure  1.  Geometry  of  microreactor  model.  All  units 
are  given  in  cm,  except  for  the  explicitly  labeled 
reactor  channel  which  is  0.5  mm  or  500  pm. 

Inside  the  reactor  channel  was  manganese 
dioxide  (Mn02),  a  known  Et202  decomposition 
catalyst6.  The  activation  energy  and  frequency 
factor,  the  kinetic  parameters  that  control  the 
reaction  rate,  were  derived  from  preliminary 
results  of  experimental  data  collected  during 
Mn02  catalyzed  hydrogen  peroxide 
decomposition.  The  activation  energy,  Ea,  and 
frequency  factor,  A,  were  20,400  J/mol  and 
2,910  mol/(m3*s*gcat)  respectively  These  values 
are  similar  to  those  in  literature  for  Mn02- 
catalyzed  H202  decomposition  reactions.7 

The  model  was  sequentially  developed  first 
solving  mass,  energy,  and  momentum  balances 
independently.  Overall  balances  were  solved 
manually  with  known  inlet  conditions  and  the 
assumption  of  full  conversion  to  verify  the 
numerical  results  from  each  of  the  independently 
solved  differential  balances.  Next,  the  model  was 
expanded  to  solve  the  differential  mass  and 
energy  balances  dependently.  Currently,  the 
model  is  expanding  further  to  include  the 
momentum  balance  to  relate  changes  in 
temperature,  concentration,  and  density. 

Initially,  the  model  treated  the  fluid  as  a 
single  phase  liquid  system.  This  preliminary 
model  lacked  an  extended  surface  and  displayed 
high  temperature  rise.  To  make  the  single  phase 
model  more  realistic,  a  stainless  steel  extended 
surface  was  added  and  the  mass  and  energy 
balances  were  refined.  Multiphase  considerations 
via  the  momentum  balance  are  gradually  being 
introduced  into  the  model.  Maintaining  the 
single  liquid  phase  assumption,  the  momentum 
was  incorporated  using  the  incompressible 
Navier  Stokes  equations.  Efforts  are  presently 
underway  to  integrate  the  non-isothermal  flow 
application  mode  to  solve  the  momentum 


balance  with  variable  density  as  a  function  of 
other  model  parameters. 

2.2  Governing  Equations  and  Assumptions 

All  equation  parameters  were  in  SI  units.  For 
specific  details  related  to  nomenclature  and  units 
of  all  constants  and  expressions  refer  to  the  table 
entitled  Model  Parameter  Details  in  the 
appendix. 

Mass  Balance.  The  mass  balance  of  the 
microreactor  system  included  the  diffusion  of 
H202  into  liquid  water  and  H202 ,  the  bulk  flow 
of  the  evolving  fluid  mixture  (consisting 
primarily  of  liquid  H202  solution  in  water,  and 
02  gas)  down  the  reactor  channel,  and 
consumption  of  H202  resulting  from  the 
decomposition  reaction.  To  model  the  mass 
balance,  the  built-in  convection  and  diffusion 
application  mode  in  COMSOL  was  used. 
Equation  1  governed  the  convection  and 
diffusion  mass  balance. 

v(-ZV/VC)+  us  VC  =  rmch  (1) 

In  equation  1,  Deff,  gave  the  diffusivity  of 
H202  into  the  liquid  water,  hydrogen  peroxide 
solution.  Due  to  time  constraints  the  diffusivity 
term  neglected  diffusion  interactions  with  the 
catalyst.  The  mass  balance  solved  for  the 
variable,  concentration,  C,  the  consumed  H202 
concentration.  The  variable  us  represented  the 
superficial  velocity  of  the  fluid.  The  rate  of 
reaction  based  on  the  consumption  of  hydrogen 
peroxide,  rH2o2 ,  was  established  assuming  first 
order  kinetics  weighted  via  an  inputted  catalyst 
mass,  Mcat  (equation  2). 

FH2O2  =  —krxnC  *  Mcat  (2) 

The  water  and  oxygen  production  rates  were 
related  to  rH2o2  (equations  3  and  4  respectively). 

ri  120  =  — ri  12O2  (3) 
r(h  =  -0.5  *rH202  (4) 

A  conventional  Arrhenius  relationship  with 
temperature,  J,  defined  the  reaction  rate 
constant,  krxm  coupling  it  to  the  energy  balance 
(equation  5).  The  activation  energy,  Ea ,  and 
frequency  factor,  A,  remained  as  previously 
defined  in  the  model  description. 

krxn  —  Ae~Ea/ RT  (5) 


Excerpt  from  the  Proceedings  of  the  COMSOL  Conference  2007,  Boston 


Energy  Balance.  The  energy  balance  of  the 
microreactor  system  included  conduction, 
advective  input  and  output  due  to  the  feed  and 
exit  streams,  exothermic  heat  generation  due  to 
reaction,  latent  heat,  and  the  heat  of  water 
vaporization.  To  model  the  energy  balance,  the 
built-in  convection  and  conduction  application 
mode  in  COMSOL  was  used.  Equation  6 
governed  the  convection  and  conduction  energy 
balance. 

V(-kefjVT)+UsCPpT  = 

(o) 

A Hrxn  *  FH2O2  -  LH  +  A Hvap  *  HfcO 

In  equation  6,  keff  gave  the  effective  thermal 
conductivity  of  the  reactor  channel  including  the 
thermal  conductivity  of  catalyst.  The  energy 
balance  solved  for  the  variable  temperature,  J, 
coupling  it  to  the  mass  balance  equation.  The 
term,  us ,  still  represented  the  superficial  velocity 
of  the  liquid.  The  term  cp  represented  heat 
capacity  and  the  term,  />,  the  density  of  the  liquid 
in  the  microchannel. 

The  heat  of  reaction,  AHrxm  defined  the  heat 
released  during  the  exothermic  decomposition 
reaction.  The  latent  heat  term,  LH ,  described  the 
contribution  of  the  heat  capacities  per  species 
required  to  increase  the  temperature  to  the 
boiling  point  of  water  (equation  7). 

LH  =  cpH2o(373  -  To)  *  rmo  +  cpo2(373  -  To)  *  ro2  (7) 

The  remaining  term  of  the  energy  balance, 
AHvap ,  expressed  the  energy  required  to  vaporize 
water  based  on  rate  of  water  production  during 
the  reaction. 

Momentum  Balance.  To  model  the 
momentum  balance,  the  built-in  incompressible 
Navier  Stokes  application  mode  in  COMSOL 
was  used.  Equation  8  (incompressible  Navier 
Stokes)  and  equation  9  (equation  of  continuity) 
governed  the  incompressible  Navier  Stokes 
momentum  balance. 

p( U  •  V)U  =  V  •  [~P2l  +  7/(VU  +  (VU)A T)]  (8) 

V  (U)  =  0  (9) 

The  density,  />,  was  modeled  as  an  average  of 
the  liquid  and  gas  comprising  the  fluid  weighted 
by  the  mole  fractions  of  H202,  H20,  and  02 
species  respectively  and  changed  as  the  reaction 
progressed  (equation  10).  This  average  density 


was  applied  under  the  assumption  that  the 
generated  oxygen  gas  was  homogenously 
dispersed  throughout  the  liquid  solution.  The 
viscosity  was  given  by  tj.  The  momentum 
balance  solved  for  the  velocity  field,  U,  and  the 
pressure,/?,  using  both  equations. 

p  -  yHiOipHiOis  +  yHiopHio  +  yoipoi  (10) 

Although,  the  simulation  results  are  not 
available  for  this  paper,  models  are  currently 
under  development  using  the  built-in  non- 
isothermal  flow  application  mode  in  COMSOL 
to  better  account  for  multiphase  effects.  The  non- 
isothermal  flow  application  solves  the 
compressible  Navier  Stokes  equations  for 
weakly  compressible  flows  (flows  with  Mach 
numbers  <  0.3). 8  Neglecting  the  effects  of  the 
catalyst  in  the  flow  channel  and  using  low  initial 
flow  rates  maintained  the  conditions  defining 
weakly  compressible  fluid  flow.  Equation  11 
(compressible  Navier  Stokes)  and  equation  12 
(equation  of  continuity  with  density  term)  govern 
the  non-isothermal  momentum  balance. 

p( U  •  V)U  -  V  •  [~p2l  +  //(VU  +  (VU)A  T)  - 
(2rj  /  3  -  /r)(V  •  U)I] 

V  •  (pU)  =  0  (12) 

The  density,  />,  remained  defined  according  to 
equation  10  previously  given.  The  momentum 
balance  solves  for  the  velocity  field,  U,  and  the 
pressure,  /?,  using  both  equations.  The  term,  1/, 
represents  the  fluid  viscosity  of  the  solution, 
whereas  k  gave  the  dilatational  viscosity.  Using 
the  solved  velocity  field,  the  mass  and  energy 
balances  couple  to  the  momentum  balance 
relating  the  density  and  viscosity  changes  to  the 
concentration  and  temperature  distributions. 

2.3  Boundary  Conditions 

For  a  tabular  synopsis  of  the  current  model’s 
boundary  conditions  and  their  affiliated 
equations  refer  to  the  table  entitled  Summary  of 
Boundary  Conditions  in  the  appendix. 

Mass  Balance.  The  only  subdomain  active 
for  the  mass  balance  was  the  reactor 
microchannel.  The  inlet  boundary  condition  was 
initial  concentration.  The  outlet  boundary 
condition  was  convective  flux.  Insulation  defined 
the  boundary  of  microchannel  walls. 
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Energy  Balance.  For  all  simulations,  the 
energy  balance’s  inlet  boundary  condition  was 
initial  temperature  and  the  outlet  boundary 
condition  was  convective  flux.  Initially,  the  wall 
boundary  was  set  to  a  constant  temperature  to 
model  active  cooling.  This  successfully  reduced 
temperature  rise,  but  remains  impractical  in 
physical  application.  We  subsequently  set  the 
boundary  condition  on  all  external  surfaces  of 
the  microreactor  system  to  heat  flux  to  simulate 
convective  cooling.  Continuity  defined  the 
remaining  interfacing  boundaries. 

Momentum  Balance.  Like,  the  mass 
balance,  the  only  active  subdomain  for  the 
momentum  balance  was  the  microchannel.  The 
boundary  condition  for  the  inlet  was  initial 
superficial  velocity.  No  slip  described  the 
boundary  condition  at  the  reactor  walls.  To 
define  the  outlet  boundary  condition  normal 
pressure/normal  flow  was  selected  and  the 
pressure  value  was  set  to  zero. 

3.  Simulation  Results 

To  obtain  a  preliminary  assessment  of 
cooling  capabilities  in  a  H202  decomposition 
microreactor,  the  microchannel  was  modeled 
alone  under  convective  cooling.  Despite 
increased  surface  to  volume  ratio,  approximately 
2500  m2/m3  compared  to  500  m2/m3  in 
conventionally  sized  reactors,  the  maximum 
simulated  temperature  was  385  K  and  the  heat 
rise  was  92  K  (Figure  2).  Even  with  the 
microchannel  geometry  under  convective 
cooling,  the  temperature  rise  was  significant. 
Figure  2  displays  the  temperature  increase, 
hydrogen  peroxide  consumption,  and  oxygen 
production  down  the  center  of  the  microchannel 
using  an  initial  volumetric  flow  rate  of  2e"9  m3/s 
of  50%  w/w  H202  for  this  simulation. 

To  facilitate  convective  cooling,  the  surface 
area  around  the  microchannel  was  extended 
Simulating  liquid  phase  mass  and  energy 
balances  in  the  microreactor  model  with  the 
extended  surface  area  reduced  heat  rise 
considerably  to  7  K.  The  maximum  temperature 
was  300  K  (Figure  3).  Figure  3  also  displays  the 
temperature  increase,  hydrogen  peroxide 
consumption,  and  oxygen  production  using  an 
initial  volumetric  flow  rate  of  2e"9  m3/s  of  50% 
w/w  H202  down  the  center  of  the  microchannel 
embedded  in  the  extended  surface  area  that  was 
convectively  cooled. 


Figure  2.  Temperature,  and  H202  and  02 
concentration  distributions  using  an  initial  volumetric 
flow  rate  of  2e'9  m3/s  of  50%  w/w  H202  down  the 
center  of  the  directly  convectively  cooled 
microchannel. 


Figure  3.  Temperature,  and  H202  and  02 
concentration  distributions  using  an  initial  volumetric 
flow  rate  of  2e'9  m3/s  of  50%  w/w  H202  down  the 
center  of  the  microchannel  including  the  convectively 
cooled  extended  surface  area. 

The  extended  surface  area  successfully  enabled 
thermal  management  with  minimal  effect  on 
oxygen  gas  generation. 

Beginning  the  integration  of  multiphase 
effects,  the  momentum  balance  was  introduced. 
Maintaining  the  single  liquid  phase  assumption 
of  the  previous  model,  the  incompressible  Navier 
Stokes  equations  were  used  to  solve  the  flow 
field  of  the  liquid  in  the  microchannel.  As 
expected  for  an  initial  volumetric  flow  rate  2e"9 
m3/s  of  50%  w/w  H202the  simulated  temperature 
rise  was  minimal  at  7  K  (Figure  4).  Figure  4 
illustrates  the  thermal  distribution  across  the 
microreactor  unit  of  this  3D  model.  The 
consumed  H202  and  02  generated  concentrations 
were  equivalent  to  the  earlier  simulations, 
indicating  modest  conversion  (Figure  5  and  6). 
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Temperature  Distribution 
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Figure  4.  Three  dimensional  temperature  distribution 
using  an  initial  volumetric  flow  rate  of  2e'9  m3/s  of 
50%  w/w  H202  over  the  convectively  cooled 
microreactor. 

Figures  5  and  6  show  the  3D  model  results  for 
the  H202  and  02  generated  concentrations 
respectively. 
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Figure  5.  Three  dimensional  H202  concentration 
distribution  using  an  initial  volumetric  flow  rate  of 
2e'9  m3/s  of  50%  w/w  H202  over  the  convectively 
cooled  microreactor. 
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Figure  6.  Three  dimensional  02  concentration 
distribution  using  an  initial  volumetric  flow  rate  of  2e" 
9  m3/s  of  50%  w/w  H202  over  the  convectively  cooled 
microreactor. 


Models  using  the  non-isothermal  application 
mode  for  the  momentum  balance  remain  under 
development,  since  the  mesh  and  solver 
parameters  require  further  adjustment  for  model 
convergence. 

4.  Summary  and  Future  Work 

Increasing  the  surface-to-volume  ratio  and 
correspondingly  the  mass  and  heat  transfer 
coefficients  via  a  microchannel  reactor  enables 
heat  management  during  catalyzed  H202 
decomposition.  To  investigate  the  mechanism  of 
thermal  management,  we  modeled  a  lone 
microchannel  reactor  (no  extended  surface)  and 
neglected  multiphase  considerations.  Significant 
heat  rise  across  the  reaction  zone  resulted  (Figure 
2).  The  surface  area  around  the  microchannel 
was  extended,  enhancing  passive  convective 
cooling  capabilities.  Although  this  secondary 
model  continued  to  neglect  multiphase 
considerations,  solving  only  the  mass  and  energy 
balances,  both  oxygen  production  and  thermal 
management  were  successfully  exhibited 
(Figures  3). 

To  integrate  the  multiphase  characteristics  of 
the  microchannel  fluid,  the  momentum  balance 
was  incorporated  and  is  gradually  being 
modified.  The  momentum  balance  was  initially 
solved  using  the  incompressible  Navier-Stokes 
equations,  which  neglected  the  flow  of  generated 
gas,  and  achieved  temperature  and  concentration 
distributions  (Figures  4-6)  that  were  equivalent 
to  the  earlier  models  (Figures  2  and  3).  The 
present  model  continues  to  demonstrate  the 
promise  of  thermally  controlled  hydrogen 
peroxide  decomposition  in  a  microreactor 
(Figures  4-6). 

The  compressible  Navier-Stokes  equations 
according  to  the  restrictions  of  the  non- 
isothermal  application  mode  in  COMSOL  are 
presently  being  incorporated  into  the  momentum 
balance.  Assuming  the  fluid  flow  is  at  low  mach 
numbers,  and  generated  oxygen  gas  is 
homogeneously  dispersed,  the  simulated 
microreactor  system  will  be  able  to  account  for 
change  in  fluid  density.  Ongoing  efforts  are 
underway  to  update  and  modify  the  model  to 
more  thoroughly  simulate  multiphase  behaviors. 
For  example,  there  are  plans  to  model  the 
influence  of  the  catalyst  on  the  diffusivity  and 
flow  field  applying  one  of  the  turbulent  flow 
application  modes.  A  comparison  between  the 
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non-isothermal  flow  application  mode  and  the 
compressible  flow  application  mode  would  also 
offer  insight  into  the  validity  of  weakly 
compressible  flow  assumption.  In  addition  to 
continuing  modeling  updates,  experimental  data 
will  be  collected  for  model  verification  and 
refinement. 

Overall,  the  modeling  results  showed  that 
thermally  manageable,  oxygen  generation  is 
viable  in  a  microreactor  provided  there  is 
adequate  external  surface  area  to  increase 
convective  cooling.  In  turn,  this  controllable 
source  of  oxygen  production  offers  a  potential 
means  for  facilitating  air- independent  SOFC 
operation. 

5.  References 

1.  UUV  Master  Plan,  Sec.  4.2.1,  pp.  60-61 
(2004) 

2.  Besser,  R.,  Microchemical  Systems  for 
UUV/AUV  Power  Systems  Presentation, 
(2004) 

3.  Papadaki,  M.;  Marques-Domingo,  E.;  Gao, 
J.;  and  Mahmud,  T.P.,  Journal  of  Loss 
Prevention  in  Process  Industrie,  06  024 
(2005) 

4.  Davis,  N.S.;  Keefe,  J.H.,  Industrial  and 
Engineering  Chemistry,  48  4  (1956) 

5.  Ehrfeld,  W.;  Volker  H.;  Holger,  L., 
Microreactors:  New  Technology  for  Modern 
Chemistry ,  1st  ed.;  pp.  5-8.  WILEY- VCH, 
Weinheim,  (2000) 

6.  Hasan,  M.A.;  Zaki,  M.  I.;  Pasupulety,  L.; 
Kumari,  K.;  Applied  Catalysis ,  181  pp. 
171-179  (1999) 

7.  Deraz,  N.  M.;  Salim,  H.H.;  El-Aal,  A.A., 
Materials  Letters ,  53  pp.  102-109  (2002) 

8.  COMSOL,  Questions  and  Answers  from  the 
Support  Desk,  COMSOL  News ,  p.17  (2007) 

6.  Acknowledgements 

The  authors  gratefully  acknowledge  the 
funding  and  resources  for  this  research  provided 
by  the  Office  of  Naval  Research  University 
Laboratory  Initiative.  We  would  also  like  to 
acknowledge  COMSOL  technical  support  for 
their  aid  throughout  the  ongoing  development  of 
this  model. 


7.  Appendices 


Table  1:  Model  Parameter  Details 

Constants 

Equation 

Symbol 

COMSOL 

assignment 

Value 

Description  (SI  units) 

vfo 

2.00E-09 

Initial  Volumetric  Flow  Rate 
(mA3/s) 

To 

293 

Initial  Temperature  (K) 

caoww 

50 

Initial  H202  %  w/w 
Concentration 

yao 

0.346 

Initial  Mole  Fraction  H202 

ybo 

0.654 

Initial  Mole  Fraction  Fl20 

rhoa 

1195 

Initial  Fluid  Density  (kg/mA3) 

P  H202s 

rha 

1097 

Average  Fl202  Density  (0- 
50%  w/w)  solution  (kg/mA3) 

P  H20 

rhb 

998 

Density  of  Water  (kg/mA3) 

P  02 

rhe 

1.3 

Average  Oxygen  Density 
between  273  and  350  K 
(kg/mA3) 

keff 

keff 

5.65 

Effective  Thermal 
Conductivity  of  Reactor 
MicroChannel  (W/m*K) 

sigmac 

3.93E-07 

Cross  Sectional  Area  of 
MicroChannel  (mA2) 

Tambient 

290 

Ambient  Temperature  (K) 

hconv 

10 

Estimated  Natural 
Convective  Fleat  Flux  for  Air 
(W/mA2*K) 

cao 

17564 

Initial  H202  Reactant 
Concentration  (mol/mA3) 

mwa 

0.034 

Molecular  Weight  of  Fl202 
(kg/mol) 

mwb 

0.018 

Molecular  Weight  of  H20 
(kg/mol) 

SHnn 

deltaH 

98200 

Heat  of  Reaction  (J/mol) 

Ea 

aE 

2.04E+04 

Activation  Energy  with  Mn02 
decomposition  catalyst 
(J/mol) 

A 

Af 

2911 

Frequency  Factor  for  Rate 
Constant 

(mol/(mA3*(s*gcat))) 

R 

Rg 

8.314 

Ideal  Gas  Constant 
(J/(mol*K)) 

radius 

0.0005 

MicroChannel  Radius  (m) 

epa 

3731 

Average  H202  (0-50%w/w) 
solution  heat 
capacity(J/kg*K) 

Mcat 

Meat 

0.03 

Mass  of  Catalyst  (g) 

Deffhh2o 

7.85E-10 

Average  effective  diffusivity 
of  H202  into  water  (mA2/s) 

Deffhh 

2.71E-12 

Average  effective  diffusivity 
of  H202  into  H202  (mA2/s) 

rhoss 

8000 

Density  of  Stainless  Steel 
316  (kg/mA3) 

cpss 

500 

Heat  Capacity  of  Stainless 
Steel  316  (J/kg*K) 

Kss 

16.3 

Thermal  Conductivity  of 
Stainless  Steel  316 
(W/(m*K)) 

CPh20 

cpb 

75.43 

Heat  Capacity  of  Water 
(J/mol*K) 

CP  02 

epe 

29.39 

Heat  Capacity  of  Oxygen 
(J/mol*K) 

Po 

1.01  E+05 

Initial  Pressure  at  Inlet  (Pa) 

dvis 

1.15E-03 

Average  dynamic  viscosity 
of  0-50%  w/w  H202  solution 
(Pa*s) 

cco 

0 

Initial  concentration  of 
oxygen  (mol/mA3) 
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Scalar  Expressions 


Equation 

Symbol 

COMSOL 

assignment 

Expression 

Description  (SI  units) 

us 

vfo/sigmac 

Initial  Superficial  Velocity  of 
Fluid  (m/s) 

k  rxn 

krxn 

Af*exp((-aE 

/(Rg*T))) 

Reaction  Rate  Constant  as 
a  Function  of  Temperature 
(mol/(mA3*s*gcat)) 

r  H202 

rt 

-krxn*Mcat*c 

Rate  Law  for  Elementary  1  st 
Order  Irreversible  Reaction 
(mol/((s*mA3)) 

r  H20 

rtb 

-rt 

Rate  Law  for  Elementary  1st 
Order  Irreversible  Reaction; 
Water  Production 
(mol/((s*mA3)) 

r02 

rtc 

Rate  Law  for  Elementary  1  st 
Order  Irreversible  Reaction; 
Oxygen  Production 
(mol/((s*mA3)) 

Xconv 

Conversion 

C  02 

co2 

cao*(0.5*Xconv) 
/(1+0.5*yao*Xco 
nv)*(T  o/T) 

Concentration  of  Oxygen 
Generated  (mol/mA3) 
assuming  negligible 
pressure  change 

Ch20 

ch2o 

cao*((ybo/yao)+ 

Xconv) 

Concentration  of  water 
(mol/mA3)  assuming 
negligible  pressure  change 

yh2o2 

yh2o2 

c/(c+co2+ch2o) 

Mole  fraction  of  H202 

yh2o 

yh2o 

co2/(c+co2+ 

ch2o) 

Mole  fraction  of  H20 

yo2 

yo2 

ch2o/(c+co2+ 

ch2o) 

Mole  fraction  of  02 

LH 

hnipc 

cpb*(373- 

To)*rtb+cpc* 

(373-To)*rtc 

Latent  Fleat  (Heat  Needed 
to  Initiate  Phase  Change) 
(W/mA3) 

A Hvap 

deltaV 

-49.004*T  + 
58754 

Heat  of  Water  Vaporization 
as  a  Function  of 
Temperature  (Heat  Needed 
to  Vaporize  Water)  (J/mol) 

Deff 

Deffhh*yh2o2+D 

effhh2o*yh2o 

Diffusion  of  H202  into 
Liquid  Solution  (mA2/s) 

P 

rhom 

c/(c+co2)*rha  + 
co2/(c+co2)*rhc 

Average  Density  of  H202 
solution  and  Oxygen  Gas 
Mixture  (Kg/mA3) 

Table  2:  Summary  of  Boundary  Conditions 

Mass  Balance 


Boundary 

Boundary  Condition 
Assigned  in 
COMSOL 

COMSOL  equation 

Inlet 

Concentration 

C=co=cao 

Outlet 

Convective  Flux 

n'(~D  effWQ  =0  t 

Reactor  Walls 

Insulation  /  Symmetry 

n  -(-D  eff^-yC+Cus)  =0 

Energy  Balance 

Inlet 

Temperature 

r=To 

Outlet 

Convective  flux 

n-(-k eir\7V  =  0 

Interfaces 

Continuity 

n(ql-q2)=0; 

qi  =-ke//;  V  Ti  + 

P  i  c  pi  Ui  Ti) 

External  Surfaces 

Heat  Flux 

ni~KffVT+PcpUT)  = 
hconv*(T -T  ambient) 

Momentum  Balance 


Inlet 

Inflow  /  Outflow 
Velocity 

U=< uso,  0,  0> 

Outlet 

Normal  Flow/Pressur^/^0 

Reactor  Walls 

|No  Slip  | 

\u=o 

f  n  represent  the  outward  normal  vector 
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STEVENS 


liHtiluk  Df  TftchnulDgv 


Crif.im  RiiHich 


Overview 

I.  Motivation  &  Objectives 


Solid  Oxide 
Fuel  Cell 


N> 


▼ 

lAnode  a  a 

Electrolyte  i 

32- 

Cathode 

h 
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II.  Initial  H202  Decomposition  Model 


III.  Current  Model 

IV.  Future  Model  Refinement 


02  feed 
from  multi-unit 
microreactor 
system 


HI 
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US  Navy’s  Unmanned  Undersea  Vehicles  (UURT 


Motivation 

H?Q?  decomposition 


2  H202 


2  H20  +  02 


Dense  Oxygen  Source 

(0.0215  moles  02/ml  in  50%  liquid  H202 

vs 

Highly  exothermic  reaction 


OITim  dF  ft  aval  Rcsnnich 


Microreactor  Features: 

4040  m2/m3  >  conventional  102  m2/m3 
improves  mass  &  heat 
transfer 


Easy  scale  up 
to 

Microchemical 
System 


o  1/5 


zzz: 


UUVs: 

Lack  efficient, 
air-independentl 
power  sources 

Fuel  cells: 

High  endurance  energy  supply;  but 

need  oxidant  in  air-independent 
environment 


Major  Modeling 
Objectives 

Display  temperature  control 

Determine  max  [H202]  for 
controlled  decomposition 
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Initial  H202  Decomposition  Reactor  Model 

Description  &  Assumptions 


CXTim  of  Naval  Rtiinmch 


[H202]  (I) 
Reactant 


Known 
Inputs  (3) 
M, 


in 


XH202in 

2CH20in 


Wall  Boundary  Condition : 
Convective  Cooling  or  Constant  Temperature 


Mn02  Catalyst 

H202(|)^H20(gi|)+1/402(gi|) 

Reaction 

^conv  =  (nH2Q2in  -  nH2Q2out) 
nH202in 


Initial  Reactor  Length,  L  (<  5cm) 

Single  liquid  phase  system  based  on 
consumption  of  liquid  H202 


Products: 
o2(g,i)  & 
H20(g,l) 


► 


Unknown 
Outputs  (fl): 


^liqout 

%H2021 


Xh201 


sout 


M 


gasout 


XmOg 

5Co2g 
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Initial  H202  Decomposition  Reactor  Model  # 

^  - —  ■■  -■-  OifiQfl  gl  H  aval  R«n  □*  ch 

Steady  State  Analysis  with  Stated  COMSOL  Application  Modes 
Governing  Equations  (all  Units  in  SI) 


Mass  Balance:  Convection  and  Diffusion 


[H202]  (variable) 

V  (— DhbCteeff  V  CH2O2)  +  Us  •  VCH2O2  =  TH202 


Average  Diffusivity 
Liquid  Solution 
(constant) 


Superficial 

Velocity 

(constant) 


1 st  order  rate  of  reaction 

rH202  =  —krxr&AiOi  *  Meat 
Rate  Constant 
(F(T)) 


Catalyst 

Mass 


Energy  Balance:  Convection  and  Conduction 


Temperature 

(variable) 


Average  Heat  of 

Density  Reaction 

(constant)  (constant) 


V(-keffVT)+  UsCppVT  =  AH-xn  *  TH202 


Effective  Thermal  Average  Heat  Capacity 

Conductivity  of  Material  of  Fluid  in  Micro¬ 
in  MicroChannel  (constant)  channel  (constant) 


E.  Lennon.  A.  A.  Burke.  R.  S.  Besser 


5 


Initial  H202  Decomposition  Reactor  Model  -# 

Simulation  FindinQS(lnitially  50%  weight  H202  concentration  for  all  runs) 


CriTiQfl  of  hiivHl  Resnmch 


00]  ftOtS  0.02  0.02s 


0  Reactor  Length  (m) 

Conclusion: 


O.CdS  0(M 


Conversion  of  reactant 
effectively  1 00%  for  reactor 
lengths  greater  than  0.04  m 
for  illustrated  volumetric 
flow  rates. 


1.8eV 


Volumetric  Flow  Rates  (m3/s): 
2e-9,  2e-10,  1e-10,  2e-11 


Conclusion: 


Heat  rise  occurs  due  to  heat 
generated  from  reaction,  even 
under  convective  cooling. 
Need  additional  cooling 
mechanism 


Reactor  Length  (m) 


0.05 
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Initial  H202  Decomposition  Reactor  Model 


Simulation  Findings:  Extending  the  Surface  Area 


M 
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H202  Concentration  (mol/m3) 


Initial  H202  Decomposition  Reactor  Model  # 

CrifiQfl  of  NilvhI  ft^srirnch 

Basis  for  Experimental  Reactor  and  Present  Model 


Note  -  reactor  is  sealed  with  a  plexiglas  cover  for  imaging  abilities 


Reactant 


Screw  holes  to  affix  plexiglas  cover 


Oring  groove  for 
cover  seal 


0.5  cm 


1  of  4  Thermo¬ 
couple  holes 


cm 


Products  Outlet 


Current  Microreactor  Model  # 

"  ~  ■"'  —  •■*  OIFiqa  gl  Naval  Rai  n  arch 

Steady  state  analysis  using  Stated  COMSOL  Application 
Modes 


Mass  Balance  (MB):  Convection  and  Diffusion 
V  (- DffeOaeffV  OY2O2 )  +  Us  •  V  Q-I2O2  =  l*H202 


Energy  Balance  (EB):  Convection  and  Conduction 

V (-keff  V  7)  +  UsCppV  T  =  AH-xn  *  mo  -  S H+  AH/ap  *  mo 


Momentum  Balance  (MoB)  3.3:  Nonisothermal  Flow 

P(U  •  V)U  =  V  •  [-p2I  +  n(VI7  +  (' V  t/)A  T)  -  (2r)/3-  k)(  V  •  t/)I] 

V  •  (pU)  =  0 
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Concentration  (mol/m3) 


Current  Microreactor  Model _ 

Intermediate  Model  Findings  (version  3.2) 


ri 


Office  at  Naval  R<«urch 


Concentration  and  Temperature  Distributions 
Down  Center  of  MicroChannel 


MicroChannel  Length  (m) 

-  -  [H202]  —  -  -  (02]  Temp 
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Temperature  (K) 


Current  Microreactor  Model _ 

Numerical  Verification  (all  Units  in  SI) 


CMTim  of  ft  aval  Rtiinmch 


Initial  Vol 

umetri 

m3 

:c  Flow  Rate  (e~9 

Zs) 

Balance 

20. 

0 

8.3 

0 

5.0 

0 

2.7 

0 

2.0 

0 

0.8 

3 

%  Nonclosure 
Energy 
Balance 

17 

11 

6 

2 

6 

18 

%  Nonclosure 
Mass  Balance 

27 

10 

6 

3 

2 

0 

%  Difference 
Reynolds  # 

17 

17 

17 

16 

15 

13 
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Current  Microreactor  Model 


ptfiui  nf  Naval  Ftasnnich 


Simiilatinn  FinHinnQ1  Tpmn  /  H _ 
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Current  Microreactor  Model  $ 

. ™  CMTSc*  gl  Naval  R^inorch 

Simulation  Findings:  X  Direction  Velocity  (m/s) 


Max:  8.128e-3 

xlO"3 
— IS.  128 

—  6.377 


5.627 


4.376 

3.126 

-1.876 

-4x625 

Min:  6.252e-4 


vfo(5)=2e-9  Isosurface:  x- velocity1  [rri/s] 
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Concentration  [mol/m  ] 


Current  Microreactor  Model 


tHFiui  of  Naval  Rcsnmch 


Simulation  Findings:  Hydrogen  Peroxide  Consumption 
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Current  Microreactor  Model 

Simulation  Findings:  Oxygen  Gas  Generation 


CHFidu  of  Naval  Rcsnnich 
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Current  Microreactor  Model 


Simulation  Findings:  Temperature  Rise 
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Simulations  and  02  Required  for  UUV  # 

L  I  ■  I  CriTic*  Dr  Naval 

Ex.  For  2.5  kW  UUV  Power  Output  Need  0.6  moles  02 
/min 

With  20  cm  tube  for  100%  conversion  of 
8.3e'10  m3/s  flow  w/  [H202]in  =  1 .8e4  mol/m3 

02  gen./tube  =  4.48e-4  mol/min 

#  tubes  for  2.5  kWsystem  =  1300 

Reactor  Volume  = 

1300  tubes  (0.16  cm3/tube)  =  200  cm3 

A0z5_Ljter  space  would  likely  be  enough  for 
this  reactor  with  cooling  components 
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Future  Model  Refinement 


CMFmfl  at  Naval  R**Mfch 


•  Improve  solver  parameters  to  minimize 
artifacts  on  concentration  profiles 

•  Continue  integrating  multiphase 
considerations 

-  Different  types  of  flow  (ie  slug  flow) 

-  Bubble  generation  from  02  gas 

•  Compare  experimental  data  to  be  collected 
for  continued  model  refinement 


IV 
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CMFmfl  □I'haval  Ra*Mfch 
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Table  1:  Model  Parameter  Details 

Constants 

Equation 

Symbol 

COMSOL 

assignment 

Value 

Description  (SI  units) 

vfo 

2.00E-09 

Initial  Volumetric  Flow  Rate 
(mA3/s) 

To 

293 

Initial  Temperature  (K) 

caoww 

50 

Initial  H202  %  w/w 
Concentration 

yao 

0.346 

Initial  Mole  Fraction  H202 

ybo 

0.654 

Initial  Mole  Fraction  H20 

rhoa 

1195 

Initial  Fluid  Density  (kg/mA3) 

P  H202s 

rha 

1097 

Average  H202  Density  (0- 
50%  w/w)  solution  (kg/mA3) 

P  H20 

rhb 

998 

Density  of  Water  (kg/mA3) 

P  02 

rhc 

1.3 

Average  Oxygen  Density 
between  273  and  350  K 
(kg/mA3) 

k, 

keff 

5.65 

Effective  Thermal 
Conductivity  of  Reactor 
MicroChannel  (W/m*K) 

sigmac 

3.93E-07 

Cross  Sectional  Area  of 
MicroChannel  (mA2) 

Tambient 

290 

Ambient  Temperature  (K) 

CMFiDB  of  Nival  Ri^Hfch 


hconv 

10 

Estimated  Natural 
Convective  Heat  Flux  for  Air 
(W/mA2*K)  | 

cao 

17564 

Initial  H202  Reactant 
Concentration  (mol/mA3) 

mwa 

0.034 

Molecular  Weight  of  H202 
(kg/mol) 

mwb 

0.018 

Molecular  Weight  of  H20 
(kg/mol)  j 

A  Hrxn 

deltaH 

98200 

Heat  of  Reaction  (J/mol) 

Ea 

aE 

2.04E+04 

Activation  Energy  with  Mn02 
decomposition  catalyst 
(J/mol) 

A 

Af 

2911 

Frequency  Factor  for  Rate 
Constant 

(mol/(mA3*(s*gcat))) 

R 

Rg 

8.314 

Ideal  Gas  Constant 
(J/(mol*K)) 

radius 

0.0005 

MicroChannel  Radius  (m) 

cpa 

3731 

Average  H202  (0-50%w/w) 
solution  heat 
capacity(J/kg*K) 

M  cat 

Meat 

0.03 

Mass  of  Catalyst  (g) 

Deffhh2o 

7.85E-10 

Average  effective  diffusivity 
of  H202  into  water  (mA2/s) 
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Deffhh 

2.71  E-12 

Average  effective  diffusivity  1 
of  H202  into  H202  (mA2/s)  | 

rhoss 

8000 

Density  of  Stainless  Steel  I 
316  (kg/mA3) 

cpss 

500 

Heat  Capacity  of  Stainless  1 
Steel  316  (J/kg*K) 

Kss 

16.3 

Thermal  Conductivity  of 
Stainless  Steel  316  f 

(W/(m*K))  I 

|<^  h20 

cpb 

75.43 

Heat  Capacity  of  Water  \ 
(J/mol*K)  j 

VP  02 

cpc 

29.39 

Heat  Capacity  of  Oxygen  1 
(J/mol*K)  'I 

Po 

1.01E+05 

Initial  Pressure  at  Inlet  (Pa) 

dvis 

1.15E-03 

Average  dynamic  viscosity  I 
of  0-50%  w/w  H202  solution  1 
(Pa*s) 

cco 

0 

Initial  concentration  of 
oxygen  (mol/mA3)  | 

OtF.Qfl  of  Naval  R<i.urch 


\  Scalar  Expressions 

Equation 

Symbol 

COMSOL 

assignment 

Expression 

Description  (SI  units) 

US 

vfo/sigmac 

Initial  Superficial  Velocity  of 
Fluid  (m/s) 

h 

rv  rxn 

krxn 

Af*exp((-aE 

/(Rg*T))) 

Reaction  Rate  Constant  as 
a  Function  of  Temperature 
(mol/(mA3*s*gcat)) 

r  H202 

rt 

-krxn*Mcat*c 

Rate  Law  for  Elementary  1st 
Order  Irreversible  Reaction 
(mol/((s*mA3)) 

r  H20 

rtb 

-rt 

Rate  Law  for  Elementary  1st 
Order  Irreversible  Reaction; 
Water  Production 
(mol/((s*mA3)) 

r  02 

rtc 

0.5*rtb 

Rate  Law  for  Elementary  1st 
Order  Irreversible  Reaction; 
Oxygen  Production 
(mol/((s*mA3)) 

Xconv 

(cao-c)/cao 

Conversion 

C  02 

co2 

cao*(0.5*Xconv) 

/(1+0.5*yao*Xco 

nv)*(To/T) 

Concentration  of  Oxygen 
Generated  (mol/mA3) 
assuming  negligible 
pressure  change 

C  H20 

ch2o 

cao*((ybo/yao)+ 

Xconv) 

Concentration  of  water 
(mol/mA3)  assuming 
negligible  pressure  change 
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\vh2o2 

yh2o2 

c/(c+co2+ch2o) 

Mole  fraction  of  H202  j 

\vh2o 

yh2o 

co2/(c+co2+ 

ch2o) 

Mole  fraction  of  H20  I 

lyo2 

yo2 

ch2o/(c+co2+ 

ch2o) 

Mole  fraction  of  02 

\lh 

hnipc 

cpb*(373- 

To)*rtb+cpc* 

(373-To)*rtc 

Latent  Heat  (Heat  Needed  1 
to  Initiate  Phase  Change) 
(W/mA3)  '1 

A  Hvap 

deltaV 

-49.004*T  + 
58754 

Heat  of  Water  Vaporization  1 
as  a  Function  of  j 

Temperature  (Heat  Needed 
to  Vaporize  Water)  (J/mol)  | 

Deffhh*yh2o2+D 

effhh2o*yh2o 

Diffusion  of  H202  into 
Liquid  Solution  (mA2/s) 

rhom 

c/(c+co2)*rha  + 
co2/(c+co2)*rhc 

Average  Density  of  H202 
solution  and  Oxygen  Gas 
Mixture  (Kg/mA3) 

Table  2:  Summary  of  Boundary  Conditions 

Mass  Balance 


Boundary 

Boundary  Condition 
Assigned  in 
COMSOL 

COMSOL  equation 

Inlet 

Concentration 

C=  co=cao 

Outlet 

Convective  Flux 

n'(-D  effXJC)  =0  t 

Reactor  Walls 

Insulation  /  Symmetry 

w(-D  eff^C+Cus)  =0 

Energy  Balance 


Inlet 

Temperature 

T= To 

Outlet 

Convective  flux 

n'(-k effSJT)  =  0 

Interfaces 

Continuity 

n(ql-q2)=0; 

qi  --Kf/i  V  Ti  + 

P  i c  Pi  Vi  Ti) 

External  Surfaces 

Heat  Flux 

ni-Kff\7T+PcpUT)  = 
hconv*(r  -T  ambient) 

Momentum  Balance 


Inlet 

Inflow  /  Outflow 
Velocity 

U=< uso,  0,  0> 

Outlet 

Normal  Flow/Pressure 

P  =Patm 

Reactor  Walls 

No  Slip 

U= 0 

f  n  represent  the  outward  normal  vector 


